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ABSTRACT 

The quasi-periodic feature of 1-10 years exhibited in pulsar timing noise has not been 
well understood since 1980. The recently demonstrated correlation between timing 
noise and variation of pulse profile motivates us to further investigate its origins. Wc 
suggest that the quasi-periodicity feature of timing noise, with rapid oscillations lying 
on lower frequency structure, comes from the geodetic precession of an unseen binary 
system, which induces additional motion of the pulsar spin axis. The resultant change 
of azimuth and latitude at which the observer's line of sight crosses the emission 
beam is responsible for the variation of timing noise and pulse profile respectively. 
The first numerical simulation to both timing noise and pulse profile variation are 
thus performed, from which the orbital periods of these pulsars are of 1-35 minutes. 
Considering the existence of the ultracompact binary white dwarf of orbital period of 
5.4 minute, HM Cancri, such orbital periods to pulsar binaries are not strange. The 
change of latitude of the magnetic moment exceeding the range of emission beam of 
a pulsar results in the intermittency, which explains the behavior of PSR B1931+24. 
Therefore, it provides not only a mechanism of the quasi-periodic feature displayed on 
some "singular" pulsars, like timing noise, variation of pulse profile and intermittency; 
but also a new approach of searching ultra-compact binaries, possibly pulsar-black hole 
binary systems. 
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1 INTRODUCTION 

Pulsar timing noise is the discrepancy between pulsar ar- 
rival times at an observatory and times predicted from 
a spin down model, containing the rotational frequency 
and its derivative. Such low-frequency structures previ- 
ously observed in p ulsar data sets have been explaine d 
by random processes (|Cordes et al.lll98"ol: Lvne et al.ll200ol ). 
unmodelled planetary compa nions ( Corded 19931 ) or free- 
precession jStairs et aLlbood ). However, these models pre- 
dict either random or exact periodic timing delay, which 
cannot well account for the quasi-periodic feature revealed, 
i.e., in the large invest igation of pulsar rotation proper- 
ties (|Hobbs et al. II2OI0I ). 

Lyne et al. (2010) reported the correlation of timing 
noise with changes in the pulse shape in six pul sars. And 
the link among mode changing (|Wang et al.l 120071 1 . nulling. 
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intermittency, pulse-shape variability, and timing noise is 
addressed. These phenomena are attributed to t he change 
in the pulsar's magnetosphere (|Lyne et al. II2OI0I ). 

This paper shows that the geodetic precession induced 
motion of pulsar spin axis can change both the latitude and 
azimuth of the magnetic moment of the pulsar, which mimics 
the magnetosphere change, and hence well account for the 
quasi-periodic structure and its correlated phenomena. 

In the absence of spin precession, the magnetic 
moment of the pulsar, n, with a misalignment angle, 

a, r elative to the spin axi s, rotates around a fixed 

axis (IRadhakrishnan fc Cook 1 119691 ). The observed pulsed 
emission occurs when the pulsar beam sweeps past the line 
of sight, in the case of dipole magnetic field, as shown in 
Fig. 1. 

If a pulsar is in orbit with a rotating companion star, 
then the two spinning tops are in the gravitation field of 
each other, and both the spin axes will precess around the 
common orbital angular momentum vector. This is the so- 
called geodetic precession. The precession of a pulsar spin 
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axis will cause additional motion of the magnetic moment, 
/X, owing to a fixed, a, the misalignment angle between the 
spin axis and ^. Consequently, both the azimuth and lati- 
tude at which the observer's line of sight crosses the beam 
change, which in turn leads to quasi-periodic timing delay 
and variation in pulse profile respectively. However, these 
pulsars with significant timing noise are usually identified 
as isolated pulsars instead of binary ones. This problem can 
be explained by three reasons. 

Firstly, the current limit on binary search is 90 min- 
utes for radio pulsars, and 10 minutes for X-ray pulsars. In 
other words, pulsars with orbital period shorter than these 
thresholds have not been detected. But from the stand of 
evolution of binary pulsars, it is difficult to understand why 
binary pulsars with orbital period less than 90 minutes or 
10 minutes cannot exist. Obviously, the limit is of technique 
rather than true existence. Therefore, it is possible that some 
of these isolated pulsars with significant timing noise are bi- 
nary pulsars beyond the threshold of binary detection. 

Secondly, e.g., PSR 1828-11 has several 8-min continu- 
ous time series observed by Parkes telescope (R.N. Manch- 
ester, private communication), but none of them cover a full 
proposed orbital period of ~20 min. Thus, orbital modula- 
tion on the timing of PSR 1828-11 has never been observed 
in over a complete orbital period, although it has been ob- 
served for years. 

Thirdly, the short-term effect of a binary system, i.e., 
the time delay corresponding to the propagation time of 
pulsar emission during its orbital motion, the Romer delay, 
depends on the amplitude of semi-major axis of a binary, 
a. Apparently a small a leads to low signal to noise ratio in 
time of arrival. Moreover, ultra-compact binaries correspond 
to rapid orbital motion which prevents from having good 
quality profile (to discern orbital modulation) in such short 
orbital bins. 

Whereas, the geodetic precession induced additional 
time delay is a long-term effect (usually lasts a few or- 
ders of magnitude longer than that of the short-term ef- 
fect), the amplitude of which can be large in the case of 
small a. As a result, the orbital effect of binary pulsars with 
small a is difficult to detect. In contrast, the corresponding 
long-term effect is still easy to find. Therefore, such pulsars 
could be identified as si ngular pulsars with strong timing 
noise (jCong II2005I . l2006h . 

How the geodetic precession invokes the quasi-periodic 
change of timing noise and pulse profile ? In the gravita- 
tional two-body system, the geodetic precession makes one 
body precess in the gravitational field of the other, which 
influences pulsar timing in two ways. 

The first one is the spin vector of one body precesses 
around the orbital angular momentum vector, the pre- 
cessi on velocity of which is gi ven by the two-body equa- 
tion (jBarker fc 0'Connelll[l975l l. 



• L(4-|-3m2/mi) ^ 2/3 
"1 = L - 



2^ 3/5 M^(4 + 3iWr) ~ 



where L and L, are magnitude and unit vector of the or- 
bital angular momentum, mi and m2 are masses of the 
two bodies, and r the separation between them; Mr = 
mi/m2, Tq = GMq/c^ = 4.925490947/xs and = Pb{l - 
e'^) 5 /(mj) 5 (where m* = mi/M©, Pt, and e are orbital 
period and eccentricity respectively). Notice that in Eq.([T|) 



the e xpression without is th e original two-body equa- 
tion (jBarker fc O'Connelf [l975h . and the one with Tq is 
rewritten for the convenience of numerical calculation. Here 
after the observed pulsar is labeled by subscript 1, and its 
companion star by 2. 

As shown in Fig.([l]), the geodetic precession can change 
the orientation of spin axis, and thus affects the azimuth and 
the latitude at which the observer's line of sight crosses the 
beam. Such an additional motion of a pulsar spin axis can 
be analyzed equivalently by assuming a constant pulsar spin 
angular momentum vector, but a varying line of sight (LOS) , 
in a reference frame defined by the spin angular momentum 
vector of the pulsar, Si, (i', j', k'). 

Another useful reference frame is defined by the angu- 
lar momentum vector, L, (i, j, k), in which the direction 
toward the observer, n, is in the plane defined by the unit 
vectors, i and k, n = (sin i, 0, cos i). Transforming from the 
reference frame defined by L, to the reference frame defined 
by Si, t he LOS is changed fr om unit vector, n to n', which 
is read l|Konacki et al. hooj ). 



n = i?(-A)i?(*)n, 



(2) 



where n first rotates around the axis k by the angle, then 
around the axis j by angle, —A. 

A pulse is identified when the magnetic moment, /x, 
crosses the plane defined by the vectors, n', and Si and, 
at the same time, the angle between n', and fj, does not ex- 
ceed the opening semi-angle p, as shown in Fig. 1. Obviously, 
the direction toward the observer n' changes with the pre- 
cession phase ^I", which can make a pulse arrive advanced or 
retarded than that of a fixed n '. This effect can be q uantified 
by the azimuthal angle of n' (|Konacki et al. |[2003l ). 



tanO(t) — n'y{t)/n'^{t) 



(3) 



Apparently, a two-body system with a single spin precessing 
as Eq.([l} will result in the timing effect of Eq.(|3]), which pre- 
dicts a periodic timing effect similar to the free precession. 

Here is the second way that the geodetic precession in- 
fluences pulsar timing. A two-body system with spin can 
cause an additional perturbation to the classical two-body 
system without spin, and hence, all six orbital elements of a 
binary vary, including uj and n, the precession of periastron 
and th e longitude of precession of the orbital plane respec- 
tively (jSmarr fc Blandfordlll976h . However, the contribution 
of UJ and Q. to Romer delay d epends on the projected semi- 
major axis, X (|Lai et al.lll995h . On the contrary, the preces- 
sion induced long-term azimuthal time delay, as shown in 
Eq.© and Eq.Q, is not directly related with x, which can 
be significant in the case of small x. 

The arrangement of this paper is as follows. Section 2 
shows the long-term additional time delay and change of 
pulse width caused by Q,, the precession speed of the orbital 
angular momentum vector, L, around the total angular mo- 
mentum, J. Section 3 applies the model to three pulsars, 
fitting their quasi-periodicity timing noise, and simulating 
the corresponding variation of pulsar profile width and in- 
termittency for the first time. The discussion and conclusion 
of results are presented in the final Section. 
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2 LONG-TERM BINARY EFFECT ON 
PULSAR TIMING 

When Cl, the precession of the orbital angular momentum 
vector, L, around the total angular momentum, J, which 
has been ignored by Eq.((2| and Eq.((3|, is considered, the 
LOS should be transformed as, 



yV{t) = 2 cos 



.1 .COS p — cos a cos(a + /?) , 



(5) 



sin a sm{a + (3) 

Notice that the width variation of Eq.©, is derived un- 
der the assumption of circular beam, the geometry of which 
needs only two parameters, half-opening angle of the beam, 
p and magnetic inclination angle, a, as shown in Fig. 1. 



n' = R{^Xi)R{^)R{XLj)Rifl)no, (4) 

where Alj is the misalignment angle between L and J, and 
no=(0, sin J, cos J) is the unit vector of LOS defined in the 
reference frame (io, jo, ko) with ko aligned with J, as shown 
in Fig. 1. Thus, putting the components n'^{t) and n'y{t) 
obtained from Eq.Q into Eq.Q, the time delay including 
to the precession of orbital plane can be obtained. 

Note that Eq.(|3| has one more frame transformation 
than that of Eq.((2]), which is realized by the rotation of 
Ri^Lj) and In Eq.©, the LOS, no, is firstly trans- 

formed to the reference frame with k aligned with L. And 
then it transforms from the frame defined by L to the 
frame defined by the spin angular momentum. Si, as that 
in Eq.©. 

The conservation of the total angular momentum, J, 
must be held during the precession of the two bodies, which 
requires that L and S (S = S i + S2) be on the oppo- 
site side of J at any mo ment (|Smarr fc Blandfor dl 1 19761 : 
iHamilton fc SarazinI flOS^ ') . In other words, the plane de- 
termined by the three vectors, J, S and L rotates around 
the vector J instantaneously. In the study of the mod- 
ulation of the gravitational wave by Spin-Orbit coupling 
effect on merging binaries, the orbital precession velocity 
satisfying the conservation of the total angular momen- 
tum, a nd the triangle c onstrai nt have been studied exten- 
sively (|Apostolatos et al...l994i : iKiddeilll995l '). 

Apparently, the precession of the two spins. Si and 
S2, can change the total spin angular momentum vector, 
S, and hence varies Xlj through the triangle formed by 
S, L an d J. The n, Xlj can simply be obtained by, Xlj = 
sm~^{^y S% + Sy/L), as shown in Fig. 1. 

Moreover, Q can be easily obtained by the precession 
phase of the two spins, = tan~^ {—Sx/Sy), where Sx = 
Si sin Ai cos rji + S2 sin A2 cos 772 and Sy = Si sin Ai sin rji + 
52sinA2 sinr;2 in which i]i = Qit+rjio, V2 ~ ^2t+ri2o (where 
?7io and rj^o are initial phases). 

Therefore, just starting from the very basic principles, 
the conservation of the total angular momentum, J, the vari- 
ation of Xlj and Q required in Eq.Q can be obtained re- 
spectively. 

Putting the precession phase of the pulsar (in the refer- 
ence frame defined by L), 5* = ryi — SI, and Aii the misalign- 
ment angle between the spin angular momentum. Si, and 
L, into Eq.©, the additional time delay can be calculated, 
which predicts a quasi-periodic variation. Notice that the po- 
lar angle of Si and S2 with respect to L is treated as equal 
to Ai and A2, the misalignment angles with J respectively, 
which are good approximations in the case of Alj ^ 1. 

Besides the azimuthal effect, the geodetic precession can 
induce latitude effect and hence vary the pulse profile as well, 
as shown in Fig.l. The impact angle between the beam axis 
and LOS is given, (3 — cos~^ n', — a, from which the pulse 
width can be obtained (|Lvne fc Manchestedll98^ '). 



3 APPLICATION TO THREE PULSARS 

Having simple modification to the geodetic precession in- 
duced timing effect, the new model is ap plied to three 
pulsars (|Lvne et al. |[201ol : iKramer et al.ll2006l i. In the fitting 
of timing residuals, we focus on whether the very basic fea- 
tures predicted by the model fit the quasi-periodic charac- 
teristics exhibited in these pulsars. 

Monte Carlo method is performed to search the best 
combination of parameters that make Eq.© fit the tim- 
ing noise of the three pulsars, the best fit parameters are 
given in Table 1. In the fitting of a timing noise structure, 
such as Fig. 2-4, 10 parameters are used, three binary pa- 
rameters, P^, Mr, I; and six parameters for the spin angu- 
lar momenta of the two bodies. Si and S2, which are their 
magnitude (51,5*2), misalignment angle with the total angu- 
lar momentum (Ai,A2), precession phases (771,772), as well as 
the parameter, c^o, representing the initial phase of 0{t) in 
Eq.(13|, as shown in Table 1. 

Although the timing noise displayed in the three pulsars 
appear different, the relatively rapid oscillations lie on lower- 
frequency structures is common for all of them, which differs 
only in relative amplitude of oscillations and in the time 
scale of quasi-periods. 

Among the three pulsars with significant timing noise, 
two pulsars, PSR B1540-06 and PSR B1828-11 have ob- 
servat ions both in timi ng residual and variation of pulse 
width (|Lvne et al. II2OI0I ) . In our model, these two effects are 
fitted by Eq.Q and Eq.lO respectively. 

As shown in Fig. ([2]) and Table 1, the fitting of both 
the timing residual and the variation of pulse width of PSR 
B 1540-06 give two solutions, the mass ratio of which are 
Mr = 2.63 and Mr = 0.51 respectively as shown in Table 1. 
This means companion mass of 3.7Mq and O.JMq respec- 
tively, in the case of pulsar mass of 7ni — 1.4Mq. As shown 
in Table 1, an orbital period of 35 minutes and compan- 
ion mass of 3.7 Mq suggests a black hole companion to the 
binary, well this solution is limited by the large projected 
semi-axis, x predicted. 

The periodicities of 250, 500 and 1000 days have been 
rep orted in the rotatio n rate and pulse shape of PSR B1828- 
11 (|Stairs et al.ll200Gl ). These periodicities have been inter- 
preted as being caused by free precession of this pulsar, even 
though it had been argued that this was not possible in the 
presence of the superfluid compo nent believed to exist inside 
neutron stars (|Lvne et al. 1120101 ). Apparently, the geodetic 
precession induced timing and profile variation avoids this 
difficulty automatically. 

The fitting of timing residual and pulse profile variation 
of PSR B1828-11 also gives two numerical solutions. One so- 
lution corresponds to precession periods of 2n/Qi = 503d 
and 2n/Q2 — 592d respectively, which accounts for the 
~500d period shown in Fig. ([3}. The discrepancy between 
them corresponds to a period of 3000d which is responsible 
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for the low- frequency structure displayed in Fig. ([3]). And the 
~250d period can be explained by the period 2n /{^i +02). 

The discrepancy between the observed and predicted 
variation of pulse width displayed in PSR B1540-06 and 
PSR B1828-11 is mainly due to the assumption of a circular 
emission beam. 

PSR B1931-I-24 behaves as an ordinary isolated radio 
pulsar for 5 to 20 days long and switches off for the next 25 
to 35 days. Such on and off repeat quasi-periodically. More- 
over, it spins down 50% faster when it is on than when it 
is off. These two unrelated phenomena can be interpreted 
by the geodetic precession too. As shown in Fig.Q, the 
timing residual of PSR B1931-f 24 can also be fitted by the 
rapid and slow oscillation of quasi-period of 7d and 45d re- 
spectively, which correspond to the precession period of the 
two bodies respectively. Although there are much fewer ob- 
servational timing residual data points than the other two 
pulsars, the correlation of timing residual with the on and 
off states, as shown in Fig. (2)), imposes strong constraints 
on the fitting parameters. The phase discrepancy between 
timing residual and the on and off states is likely due to the 
assumption of circular emission beam. 

The rapid precession of the spin axis of PSR B1931+24 
suggests that it is in a binary pulsar system with period of 
~ Iminute as shown in Table 1. Furthermore, the fitting of 
the timing residual of PSR B1931-I-24, as shown in Fig.®, 
automatically explains the spin down 50% faster when it is 
on than when it is off occurred on this pulsar. 

The quasi-period of 40d corresponds to so large a varia- 
tion in misalignment angle between the LOS and the pulsar 
spin axis that can make the pulsar emission beam outside 
the LOS. Thus explains the on-off state and hence the in- 
termittency exhibited in this pulsar. 

Comparatively, in the other two pulsars, the amplitude 
of change of j3 as shown in Fig 1, is not sufficient to move the 
emission beam out off LOS, although the geodetic precession 
varies both the timing residual and pulse profile. 



4 DISCUSSIONS AND CONCLUSIONS 

The first numerical simulation to both timing noise and 
pulse profile variation are performed, from which the or- 
bital periods of these pulsars are of 1-35 minutes. The quasi- 
periodicity feature of timing noise, with rapid oscillations 
lying on lower frequency structure, is well explained in the 
context of the geodetic precession of an unseen binary sys- 
tem, which induces additional motion of the pulsar spin axis. 

Timing noise without such a character, e.g., PSR 
B0954+08, PSR B0642-03, PSR B1822-09 and PSR 
B0950-I-08 (|Lvne et al. I I2OIOI '). cannot be fitted by our 
model. In other words, the timing noise of these pulsars may 
have other origins. 

As shown in Fig. 2-4, it clearly shows that the geodetic 
precession induced timing effect tends to be more regular, 
while the true timing noise appears more irregular. Such a 
deviation may have three reasons, (a) The timing residual 
like Fig. ([3]) is obtained by fitting the time of arrival by the 
rotation equation with parameters, v, and v (or v). The 
values of these parameters can change the shape of the tim- 
ing residual, (b) Throughout the fitting, we assume circular 
emission beam, while the actually beam can be non-circular. 



which can also lead to irregularities in both timing residual 
and pulse profile change, (c) There can be other effects like 
the quadruple moment in the companion star. Due to more 
parameters and assumptions would be introduced to account 
for such issues, we leave it in the future investigations. 

As shown in Table 1 and the pulsar mass range dis- 
cussed above, a typical companion mass of 0.25 — O.34M0, 
corresponds to the Roche lobe radius of the companion star, 
R2, 



R2 = 0.462a( 



M2 



Ml + M2 



1/3 



3 X 10^ (cm). 



(6) 



where the semi-major axis of the binary, a, is a = 1.6 x 
lO^^cm (Notice that small ap ~ lO^cm is due to the 
nearly face one orbital inclination and the small mass ra- 
tio, 1712 /m\). The Roche lobe radius of Eq.(|6} is not much 
larger than the typical radius of White dwarf (WD) star 
of 1 X 10® cm. Thus, the possibility of a WD as compan- 
ion star cannot be exclud ed. Although mo re likely candidate 
would be a strange star (|Dev et al.l 19981 ). a quark star (|Xd 
or even a black hole. 

The fitting parameters of Table 1 indicates that the 
orbital inclination angle, i ~ I, influences the amplitude of 
the predicted timing residual. Whereas, the semi-major axis 
of a binary, a, can influence the precession velocity but not 
the amplitude of predicted timing residual. For a binary with 
a precession period of a few years, the typical semi-major 
axis, a, is of ~ lO^^cm, and with a small i, the projected 
semi-major axis, x, can be a few ms. If such pulsars are 
observed in the single pulse mode at time scale of several 
or tens of the predicted orbital period as shown in Table 1, 
their ultra-compact binary nature could be revealed. 

The flrst numerical simulation to the timing residual 
and its correlated phenomena given by this paper provides 
an extremely simple scenario to the link among the tim- 
ing noise, pulse-shape variability and intermittency. The bi- 
nary assumption appears strange. However, it corresponds 
to acceptable companion stars and the absence of binary 
detection is also reasonable. The model can be tested by 
the comparison between the predicted timing residual, as 
shown in Fig. 2-4 and the observed ones in the near future. 
And the binary parameters predicted in Table 1 can also be 
directly tested by timing series of time scale longer or much 
longer than the predicted orbital periods. E.g., PSR 1828-11 
has several 8-min continuous time series observed by Parkes 
telescope, but none of them cover a full proposed orbital 
period, ~20 min. A future continuous single pulse observa- 
tion of time scale several orbital periods or a few hours is 
expected in order to test the predicted binary nature. 

This actually proposes a new way of binary searching, 
which estimates the binary parameters of a candidate pul- 
sar through fitting of its timing noise, and then test the 
predicted binary parameters via timing series of appropri- 
ate time scale. The current threshold of orbital period of 
90min of radio binary pulsars could be broken by this ap- 
proach. And pulsar black hole binary system may be found 
by this method. 
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APPENDIX A: DETAILED ROTATION 
MATRIXES 




Figure 1. Definition of the angles and three coordinate systems. 
In the case of A^j 1, which can always be satisfied in our 
calculation, the approximation, >!' = r;i — f! and Ai^ = Ai are 
adopted, where Ai^ and Ai are the misahgnmcnt angle between 
Si and L, as well as Si and J respectively. 
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Figure 2. The observed (dotted) and fitted (solid) timing resid- 
ual and pulse width variation for PSR B1540-06. The upper pan- 
els, Al and Bl, show the fitted timing residual and variation of 
pulse width corresponding to the solution with large mass ratio 
in Table 1. And the bottom panels, A2 and B2, correspond to 
that of small mass ratio in Table 1. The horizontal line indicates 
that this pulsar is always at on-state. 



The four rotation matrixes are expressed as follows 

R{n) = 



cos SI sin n \ 
— sinfi cosQ 1 ,R{\l.j) 
1/ 



/ cos * sin * \ 
= -sin* cos* ,-R(-Ai) 
V 1/ 



1 

cosAlj sin Xlj 
— sinA_Lj cosAlj 



' cos Al — sin Ai 
1 
sin Al cos Ai 



(cos f7(cos Al cos Alj sin * -I- sin Ai sin Xlj) + cos * cos Ai sin SI) 

(n'y = COS I cos * sin Xlj + sin J(cos * cos Xlj cos SI — sin * sin Q.) 
n'^ = cos J(cos Al cos Aij -I- sin * sin Al sin Alj) + sin J 

(cos S7(cos Alj sin * sin Ai — cos Ai sin Alj) + cos * sin Ai sin Q,) 

(Al) 



then expanding Eq.(|4|, the three components of n' can be 
written as: 



— cos /(— cos Alj sin Al + cos Al sin * sin Alj) + sin 7 
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6 Gong & Li 

Table 1. The fitting parameters for three pulsars 



FSR 


rms 


PJ(0.01,2) 


M,(0.01,3) 


/(0,7r) 


ci{0. 0001, 0.01) 


C2(0.001.0.5) 


1)10 (0, 2it) 


>(2n(0,2!r) 


Ai(0, it) 


A2(0,7r) 


*o(0,2ir) 


q(0, tt) 


P(0, 30°) 


B1540 00 


1.4 


0.59 


2.63 


3.10 


0.0018 


0.026 


4.25 


3.56 


1.25 


2.51 


0.0052 


1.88 


7.5° 




1.3 


0.28 


0.51 


3.12 


0.0096 


0.063 


0.78 


4.92 


0.59 


3.08 


-0.0036 


2.54 


4.3° 


B1S28 11 


13.4 


0.22 


1.21 


2.97 


0.0001 


0.08 


3.71 


3.94 


0.82 


0.94 


0.05 


2.31 


14.4° 




10.4 


0.19 


0.11 


3.09 


0.0086 


0.055 


1.39 


2.05 


0.24 


1.44 


0.036 


2.90 


7.0° 


B1931+24 


1.3 


0.016 


0.12 


3.12 


0.00010 


0.0057 


5.18 


1.48 


1.68 


0.87 


-0.016 


1.26 


10.0° 



The magnatic inclination angle a and half opening angle of emission beam, p are obtained through fitting of pulse profile variation or 
on-off states. All angels, except for p, are in radian. Pb is in hour, Mr = 'm\lm2\ and the spin angular momenta, c\ = SxjL and 
C2 = S2/L arc dimensionless. 00 is the initial phase for 0(t). The range in the parenthesis followed each parameter indicates the 

searching parameters space for the fitting process. The root mean square{rms) of fitting is given by \Jy^^'^_-,iy'i ~yj)^/^': where 
N'{= A'^ — 10 — 1) is the degree of freedom, y^. and j/j arc the fitted and observed timing residual, respectively. 




Figure 3. same as Fig.© but for PSR B1828-11. 




-4; : 
52750 52 800 52 850 52900 52950 
MJD 



Figure 4. The observed and fitted timing residual for PSR 
B1931+24. The dots are the observed timing residual and the 
solid curve represents the fitted one. The horizontal lines repre- 
sent the on states of the pulsar predicted by the model. 
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